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Abstract. We report observations of three SDSS z > 6 QSOs at 250 GHz (1.2 mm) using the 117-channel Max- 
Planck Millimeter Bolometer (MAMBO-2) array at the IRAM 30-meter telescope. J1148+5251 (z = 6.41) and 
J1048+4637 (z = 6.23) were detected with 250 GHz flux densities of 5.0±0.6mJy and 3.0±0.4mJy, respectively. 
J1630+4012 (z = 6.05) was not detected with a 3a upper limit of 1.8 mjy. Upper flux density limits from 
VLA observations at 43 GHz for J1148+5251 and J1048+4637 imply steeply rising spectra, indicative of thermal 
infrared emission from warm dust. The far-infrared luminosities are estimated to be » 10 13 Lq, and the dust 
masses » 10 8 Mq, assuming Galactic dust properties. The presence of large amounts of dust in the highest 
redshift QSOs indicates that dust formation must be rapid during the early evolution of QSO host galaxies. Dust 
absorption may hinder the escape of ionizing photons which reionize the intergalactic medium at this early epoch. 
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1. Introduction 

The search for the most distant and early galaxies has 
become a rapidly evolving field in extragalactic astron- 
omy. Optical imaging and spectroscopic surveys (Palomar 
Sky Survey, PSS; Sloan Digital Sky Survey, SDSS, York 
et al. 2000) have revealed a large number of QSOs up 
to redshifts of 6.4 (Fan et al. 2001, 2003). About 150 
high-redshift QSOs selected from these surveys were re- 
cently observed at millimeter wavelengths, detecting ther- 
mal emission from one third of them (Omont et al. 2001, 
2003; Carilli et al. 2001a), reaching out to a redshift of 5.5 
(Bertoldi & Cox 2002). Although these optically bright 
QSOs give a somewhat biased view on the relation be- 
tween the formation of stars, massive black holes, and 
galaxies in the early Universe, they complement espe- 
cially at the highest redshifts the blank field submillime- 
ter imaging surveys, which have uncovered a population 
of z > 2 dust-obscured starburst galaxies, likely to be 
spheroidal galaxies in their formation stages (Smail et al. 
1997; Hughes et al. 1998). 
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Recently, Fan et al. (2003) discovered three QSOs at 
z > 6 in the SDSS, including J1148+5251 at z = 6.41, 
the QSO with the highest known redshift. The spectra of 
these QSOs show the Gunn-Peterson quenching of contin- 
uum emission blueward of Lya, thus probing the end of 
the reionization epoch of the Universe (White et al. 2003). 
These sources provide an opportunity to study the growth 
of massive black holes and their associated stellar popula- 
tions at the end of the "dark ages" , in the earliest epochs 
of luminous cosmic structure formation. 

In this Letter, we report the detection of 250 GHz 
(1.2 mm) continuum emission from two of the SDSS z > 6 
QSOs and upper limits to their 43 GHz continuum emis- 
sion which confirm that the emission is thermal dust radia- 
tion, thus enabling us to estimate far-infrared luminosities 
and dust masses. Throughout this paper, we adopt a A- 
cosmology with H = 71 km s^ 1 Mpc^ 1 , fl\ = 0.73 and 
Q m = 0.27 (Spergel et al. 2003). 

2. Observations 

The millimeter continuum measurements were made 
in January and February 2003 using the 117-channel 
MAMBO-2 array (Kreysa et al. 1999) at the IRAM 30 
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m telescope on Pico Veleta (Spain). MAMBO-2 has a half 
power spectral bandwidth of 210 and 290 GHz with an 
effective frequency of 250 GHz. The beam size on the sky 
is 10.7 arcsec. The sources were observed with a single 
channel using the standard on-off mode with the telescope 
secondary chopping in azimuth by 32" at a rate of 2 Hz. 
For flux calibration a number of calibration sources were 
observed, resulting in an estimated absolute flux uncer- 
tainty of 15%. The total on plus off target observing time 
was 51, 128, and 68 minutes, for J1148, J1048, and J1630, 
respectively. The data were analyzed using the MOPSI 
software package. Correlated noise was subtracted from 
each channel using the weighted average signals from the 
surrounding channels. 

J1148 was imaged with MAMBO-2 using the on-the-fly 
mapping technique with chopping in azimuth by 42 arcsec. 
Skynoise was subtracted and the double beam maps were 
combined through shift-and-add. Five maps of one hour 
duration each were combined for the final image, which is 
displayed as signal/noise contours in Fig. 1. 

Continuum observations at 43 GHz of J1148+5251 and 
J1048+4637 were done using the VLA in the D config- 
uration (max. baseline = 1 km). The sources were ob- 
served for 4 hours and 0.7 hour, respectively. Standard 
amplitude calibration was performed using 3C 286. Fast 
switching phase calibration was employed using celestial 
calibrators within 3° of the target sources. The calibration 
cycle time was 200 seconds, and the phase stability was 
excellent, with typical changes in antenna-based phase so- 
lutions between calibration scans < 10° on the longest 
baselines. Images were generated using the deconvolution 
task IMAGR in AIPS, and CLEANed to residuals of 1.5a. 
The Gaussian restoring CLEAN beam was ~ 1" (FWHM). 
The r.m.s. noise level in the final image for J1048+4637 is 
0.37 mJy/beam, and 0.11 mJy/beam for J1148+5251. 

3. Results 

The results from the MAMBO and VLA observations are 
summarized in Table □ Two of the QSOs, J1148+5251 
and J1048+4737, are detected above 5a significance at 
250 GHz (rest frame 165 fim). These QSOs are not de- 
tected at 43 GHz (rest frame 970 fim) , with upper limit 
flux densities below 1 mJy. The steep rest-frame submil- 
limeter spectral index suggests that the emission is ther- 
mal dust radiation and not synchrotron radiation. 
J1148+5251. The highest redshift source yet discov- 
ered (z = 6.41), J1148+5251 is an extremely luminous 
QSO powered by a massive black hole of 3 x 10 9 M Q ac- 
creting close to its Eddington limit (Willott et al. 2003). In 
a if'-band Keck image (Fan et al. 2003), the source (K' 
= 16.9) is unresolved and no other optical source with 
K' > 21 is present within 10" of the QSO, thereby rul- 
ing out strong gravitational lensing on arcsecond scales. 
Jl 148+5251 is detected at 250 GHz in the pointed mea- 
surements at 4.8 ± 0.8 mJy, and in the map at 5.3 ± 1.0 
mJy. It was not detected at 43 GHz with the VLA to a 3a 
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Fig. 1. The 6 arcminute field surrounding the QSO 
SDSS J1148+5251. Coordinates are offsets from the opti- 
cal QSO, which is visible as a faint blue dot. Color image: 
SDSS z band image, smoothed to 5" and logarithmically 
scaled. Red-white contours: MAMBO-2 250 GHz signal to 
noise map smoothed to 13". The contours correspond to 
2 and 4 a. The r.m.s. noise level, cr, in the proper map is 
0.9 mJy in the central 100" and rises to ~ 1.7 mJy at a 
radius 200". White contours: VLA NVSS 1.4 GHz image. 
The beam size is 45", and contour values are 2, 4, 8, 32, 
64 mJy/beam. 



limit of 0.33 mJy, implying a lower limit to the spectral 
index between 43 and 250 GHz of +1.6. 

J1148+5251 was not detected at 1.4 GHz in the VLA 
FIRST survey (Becker et al. 1995). However, the survey 
shows two bright radio sources within 1' of J1148+5251 
(Fig. 1 and Table 1). The one located north-east of the 
QSO, FIRST-NE, is also detected (unresolved) at 43 GHz, 
implying a falling spectrum between 1.4 and 43 GHz of 
index -1.2. The other FIRST source (FIRST-SW) is lo- 
cated south-west of J1148+5251, coincident with a large 
elliptical galaxy at redshift 0.05 (Fan et al. 2003), and 
is not detected at 43 GHz. The lower resolution VLA 
1.4 GHz NVSS and Westerbork 327 MHz Northern Sky 
Survey (WENSS) show that FIRST-SW is located at the 
center of extended radio emission reaching about 4 ar- 
cminutes SE-NW (Fig. 1). Note that for a 1' field one 
expects only 0.003 sources by chance with £1.4 > 8 mJy 
(Fomalont et al. 2003), so the presence of several such 
objects in the vicinity of the QSO is remarkable. 

The MAMBO image of Jl 148+5251 reveals at least 
two other sources (Fig. 0: the z — 0.05 elliptical galaxy 
(FIRST-SW) south-west of J1148+5251, and a source to- 
ward the cast (mm-E) , with no optical counterpart in the 
SDSS images. We also notice a 2 a peak which coincides 
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Table 1. Properties of the observed QSOs and of field sources near J1148+5251 
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NOTE - The optical properties are from Fan et al. (2003) and the 1.4 GHz data from the VLA FIRST survey. Upper limits are 
given at the 3<r level. * Average of on-off and map measurements. " Source may be over-resolved in the FIRST survey. 
To estimate the far-infrared luminosities, Lfir, we adopt a dust temperature of 45 K and an emissivity index /3 = 1.5, which 
is typical for the spectrum of an infrared- luminous galaxy. Increasing (3 to 2 would raise the luminosities by ~ 20%, whereas 
varying the dust temperature by ±10 K would change the luminosities by about a factor 2 down or up, respectively. 



with the strong FIRST source north-east of J 1148+5251, 
corresponding to the fainter part of a double compact op- 
tical galaxy. Several potential millimeter sources at the 3 
a level are found in the 1.2 mm map, but considering the 
size of the map they are not very significant. 

It is peculiar that the QSO is surrounded by two 
strong millimeter sources. The MAMBO deep field surveys 
(Bertoldi et al. 2000a,b; Carilli et al. 2001b) show an av- 
erage surface density of sources with flux density > 4 mjy 
of 0.02 arcmin~ 2 . The probability to find two such mil- 
limeter sources within a 2 arcmin radius from the QSO is 
only 6%. The optical images show that the QSO falls into 
a region with an overdensity of foreground galaxies sur- 
rounding the z = 0.05 elliptical galaxy. That the central 
cD galaxy of a cluster can show noticeable submillimeter 
emission was pointed out by, e.g., Edge et al. (1999). The 
association of the QSO with a millimeter-bright cD could 
thus be interpreted as due in part to a lens amplification 
of the QSO by the cluster or the dark matter halo associ- 
ated with the elliptical. However, the low redshift of the 
cD or cluster would not produce a strong amplification. 
The presence of an intervening Civ absorption system at 
z = 4.95 (White et al. 2003) hints at the possible ex- 
istence of another possible lens, but a high-redshift lens 
would not produce a large amplification either. 

It remains unclear whether the eastern millimcter 
source mm-E or the strong radio source FIRST-NE are 
associated with the QSO or with the foreground cluster. 

If we assumed that the source mm-E is lens amplified 
by a factor 2, then the chance probability to find such a 
source within 2' of the QSO is of order unity. Although 
the statistical evidence is weak, the optical, millimeter, 
and radio data hint at a mild lens amplification toward 
J1148+5251, which would lower its implied luminosity and 
dust mass. Alternatively, they hint at a possible overden- 
sity of objects near the z = 6.4 QSO. 



J1048+4637. At z = 6.23, this optically very luminous 
BAL QSO is the third most distant quasar identified to 
date. There is no evidence for arcsecond scale gravitational 
lensing of this source (Fan et al. 2003). J1048+4637 is 
detected at 250 GHz, but not at 43 GHz (Table 1), and 
we find a lower limit to the spectral index between 43 and 
250 GHz of +0.6. This QSO was not detected at 1.4 GHz in 
the FIRST survey to a 3(7 upper limit of 0.43 mjy beam -1 . 
There is, however, a possible detection of a faint, 0.46 mjy 
radio source 20" west of the QSO. 
J1630+4012. This z = 6.05 QSO is the optically 
faintest z > 5.7 QSO found in the SDSS (Fan et al. 
2003). J1630+4012 is neither detected at 250 GHz nor 
at 1.4 GHz, and no radio sources are found within 1' from 
the radio source to this limit. 

4. Discussion 

All five QSOs known to date at z > 6 (Fan et al. 2001, 
2003) were observed at 250 GHz to r.m.s sensitivities ~ 1 
mjy (Petric et al. [2003] place upper limits for J1030+0524 
and J1306+0356 at z = 6.28 and 5.99, respectively). The 
detection of two quasars reported in this Letter is consis- 
tent with the 30% detection fraction of QSOs in the red- 
shift range z = 2 to 6 surveyed at 250 GHz to mjy sensi- 
tivities (Omont et al. 2001, 2003; Carilli et al. 2001a). The 
fraction of optically luminous QSOs that are also infrared 
luminous is therefore roughly constant with redshift, out 
to the highest rcdshifts explored. If the dominant dust 
heating mechanism is radiation from young stars, the im- 
plied star formation rates in J1048+4637 and J1148+5251 
are 2000 M Q yr" 1 and 3000 M Q yr" 1 , respectively, compa- 
rable to what was derived for the 1.5 < z < 5.5 QSOs 
detected at millimeter wavelengths. 

In the optical spectra of the two z > 6 QSOs detected 
at 250 GHz, the Lya+[Nv] emission lines are relatively 
weak, in contrast to the three non-detected QSOs, which 
show stronger and sharper lines (Fan et al. 2001, 2003). 
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This trend agrees with the results of Omont et al. (1996) 
that luminous high redshift QSOs with weak broad emis- 
sion or broad absorption optical lines tend to have stronger 
millimeter emission. 

From the infrared luminosities, we derive dust masses 
of 4xl0 8 and 7x10 s M Q for J1048+4637 and Jl 148+5251, 
respectively. Following Omont et al. (2001) we here 
adopted a dust absorption coefficient at 230 ^m of K230 = 
7.5 cm 2 g" 1 , a value that applies to a galactic dust compo- 
sition and is unknown for high redshift sources. The dust 
mass is also affected by the exact temperature of the warm 
dust component and by the possible presence of an addi- 
tional cold dust component (see the discussion in Omont 
et al. 2001). Despite these large uncertainties, the esti- 
mated dust masses are huge, implying a high abundance of 
heavy elements at z « 6. This is consistent with the super- 
solar metallicities found in the three QSOs discussed here 
(Fan et al. in preparation), and with the Fe/Mg abun- 
dance ratios near or above the solar value measured in 
three other QSOs at 5.7 < z < 6.3 (Freudling et al. 2003). 

The presence of large amounts of dust at redshift 6.4 
implies that efficient dust formation took place between 
the corresponding cosmic time and the epoch of early 
reionization [z r « 17, Kogut et al. 2003), a time span 
of « 0.7 Gyr. At a constant formation rate this implies a 
net dust production rate of « 1 M Q yr _1 . 

A time span of 0.7 Gyr is short by at least a factor 2 to 
efficiently produce refractory grains in the quiescent winds 
of low- mass (M < 8M Q ) stars. If the observed dust were 
the product of stellar processes, the initial refractory dust 
enrichment might have occurred primarily through dust 
condensation in supernova remnants, and perhaps in the 
winds of high-mass (M > 40 M ) stars, which are thought 
to have dominated the early phases of star formation (e.g., 
Bromm & Loeb 2003). The dust in the early Universe 
must then be composed of silicates and perhaps oxides, 
since carbon dust is primarily formed from stars of mass 
2-5 Mq (Dwek 1998) - except if dust production in the 
winds of high-mass stars was important (Todini & Ferrara 
2001). For silicate and iron oxide dust the mass absorption 
coefficient may take a higher value (Hcnning & Muschkc 
1997) than that we adopted to compute the dust mass, 
which may therefore be overestimated. 

Elvis et al. (2002) proposed that dust may be produced 
in an outflow from the broad line region of the AGN. The 
dust content of QSOs should then be roughly independent 
of redshift, which is consistent with the similar FIR prop- 
erties of QSOs from redshift 2 to beyond 6. A difficulty 
with this model is that it requires pre-existing heavy ele- 
ments in the interstellar medium, and hence prior star for- 
mation may be required regardless. This mechanism can 
produce up to 10 7 M Q of dust fairly readily, but the pro- 
duction of the much larger amounts may be problematic. 

Although the quasars we observed are extreme and 
rare objects hardly representative of the dominant star 
forming galaxies in the early Universe, they do show that 
early star formation lead to a rapid metal and dust enrich- 
ment of the interstellar medium. Absorption by dust could 



have significantly reduced the escape fraction of ionizing 
radiation from galaxies during the epoch of reionization. 
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